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Recently,  there  has  been  considerable  interest  in  atoms  confined  to  optical  cavities  or  embedded  within 
photonic  band-gap  materials,  where  the  spectrum  of  vacuum  fluctuations  is  modified.  As  a  first  step  in  extend¬ 
ing  traditional  semiclassical  methods  into  this  regime,  the  present  paper  describes  a  semiclassical,  density- 
matrix  methodology  for  treating  radiative  processes  in  a  low-2  cavity.  This  semiclassical  approach  is  tested  by 
comparing  its  predictions  to  the  optical  cavity  experiments  of  Heinzen  and  Feld  [Phys.  Rev.  Lett,  59,  2623 
(1987)]  and  by  calculating  transient-nutation  decay  in  free  space  induced  by  a  broadband  laser.  The  method¬ 
ology  is  then  used  to  examine  transient  nutation  in  a  low-2  optical  cavity,  where  it  is  found  that  the  detuning 
between  the  atomic  and  cavity  resonances  has  a  significant  effect  on  the  transient  nutation's  decay  rate. 
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L  INTRODUCTION 

Although  the  semiclassical  description  of  radiative  pro¬ 
cesses  is  only  a  surrogate  for  quantum  electrodynamics, 
semiclassical  methods  nonetheless  remain  a  vital  tool  for 
studying  field/atom  interactions  [1].  In  particular,  processes 
ranging  from  magnetic  resonance  [2]  to  optical  pumping  [3] 
to  multiphoton  ionization  [4]  are  all  well  described  using  the 
semiclassical  density  matrix.  Notwithstanding  their  present 
utility,  however,  semiclassical  methods  typically  limit  them¬ 
selves  to  situations  in  which  a  simple  spontaneous-emission- 
decay  term  may  be  added  to  the  appropriate  equations  in  the 
fashion  of  a  Weisskopf-Wigner  approximation  [5],  While  this 
procedure  is  valid  for  treating  spontaneous  emission  in  free 
space,  since  the  correlation  time  of  vacuum  fluctuations  in 
that  case  is  infinitely  short;  it  is  at  best  problematic,  and 
generally  just  wrong,  when  considering  radiative  processes 
in  cavities  or  photonic  bandgap  materials  where  the  atom's 
vacuum  environment  is  modified. 

With  the  growing  interest  in  radiative  processes  taking 
place  within  structured  vacuum  environments  [6],  it  seems 
only  natural  to  attempt  an  extension  of  semiclassical  meth¬ 
ods  into  this  regime.  For  example,  given  recent  demonstra¬ 
tions  of  lasing  in  photonic  band  gaps  [7],  there  could  be 
value  in  describing  photonic  band-gap  lasers  by  using  stan¬ 
dard,  semiclassical  laser  theory  [8].  Additionally,  with  an  eye 
on  atomic  clock  applications,  fully  semiclassical  methods 
could  be  used  advantageously  to  describe  the  dynamics  of 
atoms  confined  in  high-2  optical  cavities.  Towards  the  even¬ 
tual  goal  of  developing  these  computational  techniques,  the 
present  work  incorporates  a  semiclassical  treatment  of  spon¬ 
taneous  emission  into  the  density-matrix  equations  in  order 
to  obtain  a  completely  semiclassical  methodology  for  de¬ 
scribing  radiative  processes  in  a  colored  vacuum  (i.e.,  a 
low-2  cavity).  The  principal  question  to  be  addressed  here  is 
whether  or  not  such  a  semiclassical  approach  is  viable  given 
the  serious  problems  encountered  in  previous  attempts  [9]. 

The  difficulty  with  the  present  program  is  that  spontane¬ 
ous  emission  is  fundamentally  different  from  other  relaxation 
processes  [10],  since  transitions  are  only  allowed  from 
higher  to  lower  energy  eigenstates.  This  “microscopic  re- 
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versibility  issue"  is  a  major  obstacle  to  ail  semiclassical 
theories  of  spontaneous  emission  [11],  and  in  the  present 
work  we  circumvent  the  problem  by  attributing  spontaneous 
decay  to  two  complementary  Hermitian  processes.  For  the 
first  of  these  processes,  the  vacuum  of  QED  is  approximated 
by  a  classical  zero-point  field  (zpf)  [12],  which  induces  di¬ 
pole  transitions  between  atomic  states.  Additionally,  so  as  to 
retain  consistency  with  QED  [13],  we  assume  that  the  atom's 
interaction  with  the  zpf  causes  it  to  emit  a  radiation-reaction 
field,  which  also  induces  dipole  transitions.  For  higher- 
energy  to  lower-energy  atomic  transitions  these  two  path¬ 
ways  add  constructively  to  produce  spontaneous  decay,  while 
in  the  case  of  lower-  to  higher-energy  transitions  the  two 
pathways  destructively  interfere,  eliminating  any  possibility 
of  “spontaneous  excitation"  [14], 

In  what  follows,  the  density-matrix  equations  are  obtained 
by  averaging  over  the  zero-point  field  and  by  assuming  that 
the  atom  does  not  modify  the  zpf  spectrum.  This  latter  as¬ 
sumption  restricts  the  present  theory  to  situations  where  the 
correlation  time  associated  with  vacuum  fluctuations  is 
shorter  than  the  atom's  excited-state  lifetime.  The  resulting 
density-matrix  equations  are  therefore  reminiscent  of  a 
master-equation-type  approach  to  the  problem,  where  the 
equations  are  obtained  by  tracing  over  reservoir  states  [15], 
However,  the  Markov  approximation  is  not  invoked,  so  that 
in  the  case  of  strong  fields  the  Rabi  period  can  be  equal  to  or 
smaller  than  the  zpf  correlation  time. 

Following  a  derivation  of  the  density-matrix  equations  in 
Sec.  II,  the  semiclassical  approach  is  tested  in  Sec.  Ill  by 
comparing  its  predictions  of  spontaneous  emission  and  the 
Lamb  shift  in  an  optical  cavity  to  the  measurements  of 
Heinzen  and  Feld  [16].  The  strong-field  manifestation  of  the 
theory  is  then  examined  by  computing  an  atomic  popula¬ 
tion's  transient-nutation  decay  in  free  space.  In  order  to  make 
this  test  a  bit  more  difficult,  we  consider  a  situation  in  which 
the  nutation  is  induced  by  a  stochastic  laser  field  (i.e.,  a 
phase-diffusion  field),  and  then  compare  our  results  to  well- 
established  expectations.  Finally,  in  Sec.  IV  transient  nuta¬ 
tion  in  a  colored  vacuum  is  examined,  and  it  is  shown  that 
the  transient-nutation  decay  rate  is  strongly  dependent  on  the 
detuning  between  the  cavity  and  atomic  resonances.  While 
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this  latter  result  is  new,  it  is  nonetheless  consistent  with  the 
findings  of  Keitel  et  al.  [17],  who  employed  QED  to  examine 
resonance  fluorescence  in  a  colored  vacuum. 

II.  DENSITY-MATRIX  EQUATIONS 

A.  Seniiclassical  spontaneous  emission 

For  simplicity,  we  consider  a  two-level  atom,  where  |l)  is 
the  ground  state  and  |2)  the  excited  state,  and  the  energy 
separation  between  the  two  levels  is  h<x>2\  •  The  perturbation 
acting  on  the  atom  is  made  up  of  three  Hemutian  terms: 
the  perturbation  due  to  dipole  coupling  between  the  atom  and 
some  real  (e.g.,  laser)  field;  the  perturbation  due  to  di¬ 
pole  coupling  between  the  atom  and  the  classical  zero-point- 
field  constmct;  and  the  perturbation  arising  from  the 

atom's  interaction  with  its  own  radiation-reaction  (RR)  field 
[18].  Thus,  we  write  the  Liouville  equation  for  the  atom  s 
density  matrix  in  the  interaction  picture  cr  as 

0-=  -  ^[(  v^+  v^+vy^).a]  =  &^+  0-®+  (1) 

where  and  Hq  is  the  unperturbed 

atomic  Hamiltonian.  We  note  that  by  associating  the  vacuum 
and  radiation-reaction  fields  with  Hermitian  operators  in  Eq. 
(1),  we  are  in  effect  following  a  prescription  of  Dalibard, 
Dupont-Roc,  and  Cohen-Tannoudji  [19].  These  authors  sug¬ 
gest  that  even  though  any  ordering  of  creation/annihilation 
operators  is  valid  in  QED,  a  symmetric  ordering  provides  the 
greatest  physical  meaning,  since  vacuum  fluctuation  and 
radiation-reaction  terms  in  QED  are  then  separately  Herrmt- 
ian.  Such  an  ordering  attributes  one  half  of  an  atom  s  spon¬ 
taneous  decay  rate  to  interaction  with  the  vacuum  and  the 
other  half  to  interaction  with  the  radiation-reaction  field, 
similar  to  what  is  indicated  by  Eq.  (1). 

Density-matrix  evolution  due  to  the  real  field  is,  of 
course,  well  known  and  leads  in  the  rotating-wave  approxi¬ 
mation  to  the  standard  radiative  equations  without  relaxation 
terms: 


(2a) 

id 

(2b) 

Here,  and  in  what  follows,  we  consider  the  real  field  as  aris¬ 
ing  from  a  linearly  polarized  laser,  so  that  is  the  detuning 
between  the  laser  and  the  atomic  resonance  (i.e.,  ^ 

is  the  Rabi  frequency  associated  with  the  laser/atom  interac¬ 
tion,  and  the  diagonal  components  of  the  density  matrix  have 
been  normalized.  Note  that  the  atom  s  interaction  with  the 
real  field  defines  a  quantization  axis  for  the  atoms  that  we 
take  as  the  x  axis  of  our  coordinate  system. 

Turning  to  the  atom's  interaction  with  the  zero-point  field, 

the  zpf  electric  field  component  E^\  is  written  as  a  sum  over 
modes  and  polarizations  X  [14]: 


Trfi 


1/2 


sK 


X[ZsK  exp(  -  / (o,t)  -  zfx  exp( i ]. 


(3) 


where  L  is  some  large  volume  (eventually  representing  free 
space  and  not  the  volume  of  some  experimental  optical  cav¬ 
ity),  is  a  polarization  vector,  is  the  frequency  of  the 
sth  mode,  and  the  z.vx  are  independent,  mean  zero,  complex 
random  variables: 

-  [14].  C,  is  a  relative  modal  amplitude  for  the  zpf 

associated  with  its  coloring  by  an  optical  cavity  as  will  be 
discussed  more  fully  below. 

As  mentioned  above,  the  atomic  dipole  moment  jl  is  ori¬ 
ented  in  space  due  to  the  atom  s  interaction  with  the  laser 
field.  Thus,  the  zpf  perturbation  on  the  atom  becomes 


1/2 


sk 


X  \/^['.v\  exp(/ (4) 

In  standard  fashion,  Eq.  (4)  may  be  employed  in  the  Liou¬ 
ville  equation  in  order  to  obtain 

/  TT 

o'22=-2|a‘2i|(^]  2 


X  Re{(Ti22;jx  exp(  -  /  <5/)}  (5) 

for  the  evolution  of  the  excited-state  density-matrix  element, 
where  o);)),  and 

/  \  1/2 


Xexp(/ S,^t)(a[\-<J22) 


(6) 


for  the  coherence  term  of  the  density  matrix.  Note  that  in 
writing  Eqs.  (5)  and  (6),  rapidly  rotating  terms  have  been 
ignored  as  they  effectively  average  to  zero  over  the  time 
scale  for  density-matrix  evolution. 

Proceeding  to  the  radiation-reaction  field's  contribution  to 
the  evolution  of  the  density  matrix,  we  have  — 

-  E^.  Here,  is  the  radiation-reaction  field  emitted  by  the 
atom,  which  almost  by  definition  is  an  atomic  operator.  In 
Ref.  [14],  the  radiation-reaction  field  was  written  as  an  Her¬ 
mitian  operator, 

(7) 


Though  this  particular  expression  for  E^  leads  to  a  non- 
Hemiitian  perturbation,  so  that  some  modified  form  would 
be  required  in  the  present  analysis,  it  is  important  to  note  that 
the  radiation-reaction  field  operator  of  Eq.  (7)  is  not  gauge 
invariant  [20].  This  may  be  immediately  recognized  through 
the  field's  dependence  on  [21].  Thus,  even  though  each  of 
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the  terms  in  Eq.  (1)  has  a  clear  physical  interpretation  (since 
they  are  all  Hermitian),  this  interpretation  would  seem  to 
depend  on  the  choice  of  gauge. 

Though  one  might  be  tempted  to  view  the  gauge  depen¬ 
dence  of  the  radiation-reaction  field  as  problematic  for  the 
present  semiclassical  approach,  in  actuality  it  should  be  seen 
as  a  requirement  for  its  consistency  with  QED.  Since  vacuum 
fluctuations  and  radiation  reaction  are  really  just  “two  sides 
of  the  same  coin"  in  quantum  electrodynamics  [13],  any  ap¬ 
pearance  of  separability  between  the  vacuum  and  radiation- 
reaction  fields  must  be  illusory.  In  the  present  semiclassical 
approach,  this  manifests  itself  as  a  gauge  dependence  of 
Thus,  there  is  no  radiation-reaction  field  that  exists  in¬ 
dependently  of  the  vacuum  in  the  semiclassical  approach;  it 
is  only  the  present  choice  of  gauge  that  allows  for  such  an 
interpretation. 

Since  (and  hence  V^)  is  not  gauge  invariant,  it  is  not 
a  measurable  quantity;  it  is  simply  a  computational  tool  [20]. 
Thus,  while  we  could  attempt  to  devise  a  colored  vacuum 
analog  to  Eq.  (7),  and  then  compute  the  perturbation's  matrix 
elements  as  was  done  above  for  the  zpf,  there  is  really  no 
point  using  that  procedure  if  a  more  expedient  approach  may 
be  found.  Here,  we  take  advantage  of  the  fact  that  on  aver¬ 
age  radiation  reaction  and  the  zpf  must  constructively  inter¬ 
fere  to  produce  excited-state  decay,  while  they  must  destmc- 
tively  interfere  to  inhibit  spontaneous  excitation  of  the 
ground  state  [14,22].  Specifically,  due  to  the  linearity  of 
density-matrix  evolution,  we  assume  that  we  can  write  the 
zpf-averaged  evolution  of  any  density-matrix  element  as 

{&%U{={f(oru))+{sicm)),  (8) 

where  /  and  g  are  arbitrary  functions  or  linear  operators. 
Then,  by  appealing  to  the  fluctuation-dissipation  theorem 

[23] ,  we  make  the  ansatz  that 

so  that 

(10) 

Evidence  for  Eq.  (10)  derives  from  the  fact  that  it  has  the 
proper  form  to  describe  spontaneous  decay  in  free  space  [14] 
and  that  it  is  consistent  with  the  QED-derived  relationship 
between  radiation  reaction  and  the  electromagnetic  vacuum 

[24] .  In  the  present  work,  we  obtain  further  validation  for  Eq. 
(10)  a  posteriori^  when  we  compare  the  predictions  of  the 
semiclassical  approach  with  experiment  in  Sec.  III.  (Further 
discussion  of  this  ansatz  is  provided  in  the  Appendix.) 

Proceeding,  Eqs.  (5)  and  (6)  are  averaged  over  the 
yielding 

/  TT  \ 

(‘^22)zpr“  j  2  ^ ^ sk) 

X  Re{{ori2<':,v\)zpi  exp(  ~  i <?/)}  (11a) 

and 


X(;;';(«^ll”Cr22)>zpf-  (11b) 

(In  what  follows  we  will  drop  the  zpf  subscript  on  the  angu¬ 
lar  brackets,  it  being  understood.)  Considering  for  the  mo¬ 
ment  just  the  coherence  term,  we  immediately  obtain  from 
Eq.  (10) 

x(:fxO-22)-  (12) 

In  order  to  evaluate  the  correlation  between  the  zpf  and 
the  density  matrix  in  Eq.  (12),  we  recognize  from  Eq.  (1)  that 

(13) 

In  Eq.  (13),  we  first  note  that  the  correlation  between  any 
single  zpf  mode  and  the  density  matrix  should  be  small, 
since  density-matrix  evolution  is  influenced  by  a  large  num¬ 
ber  of  independent  zpf  modes.  Thus,  we  make  a  (second- 
order)  decorrelation  approximation  on  the  right-hand  side  of 
Eq.  (13),  such  that  for  some  function  of  the  zpf  random  am¬ 
plitudes  '^V{Zsk)  under  the  integral  sign  of  Eq.  (13) 

(^P(:.x)«-«^)=<'I^(^u))(‘^a^)-  (14) 

We  then  get  from  Eqs.  (2)  and  (14)  In  this 

same  spirit,  we  invoke  a  mean-field  approximation  with  re¬ 
gard  to  the  atom's  radiation-reaction  field  [25],  basically  as¬ 
suming  that  due  to  the  large  number  of  zpf 

modes  that  contribute  to  this  field.  In  combination  with  the 
decorrelation  approximation,  this  mean-field  approximation 
implies  that  so  that  In  this  way, 

the  (first-order)  correlation  between  the  density-matrix  ele¬ 
ments  and  the  zpf  is  given  by  the  atom's  direct  coupling  to 
the  zpf.  Thus, 

(15) 


so  that  Eq.  (12)  becomes 


Using  Eq.  (5),  we  are  finally  led  to 
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<  o-?,) + ( o-f^) = -  ^  1 M21 1  ‘  S  \ 


(17) 

Employing  a  similar  line  of  argument  for  Eq.  (I la),  we 
obtain 

I  ^  \in 

(  ^22)  “  “  I  ^2  ]  i  (  2^3  I  2  *  ^jx) 


X  Rej  exp(  -  /  S,r )  |  ^ '  1 1 . 

which  with  the  aid  of  Eq.  (6)  becomes 


(18) 


.vX 


f  [(o-,|(/'))-(o"i->(f'))]cos[<5.v(/-r')]<3'r'. 
Jo 


(19) 


Then,  employing  Eq.  (10),  we  obtain 


477 


( d-",)  +  =  -  777 1 M21 1 '2  C; (x  •  £, J- W, 

IlL  sK 

X  fV22(''))cos[<5,,(r-/')]^J/'.  (20) 


X 


FIG.  1.  Illustration  showing  cavity  geometry  and  coordinate 
axes. 

ensemble  of  atoms  placed  in  the  center  of  a  confocal  Fabry- 
Perot  etalon.  Remembering  that  in  the  present  approximation 
to  QED,  the  zpf  is  a  real  field,  it  is  fairly  straightforward  to 
describe  the  effect  of  the  cavity  on  the  zpf:  as  the  cavity  is 
assembled,  some  of  the  zpf  modes  of  free  space  suffer  con- 
stnjctive  and  destructive  interference  as  they  reflect  off  the 
mirror  surfaces.  For  those  modes  with  |cos(^j|^cos(^c)? 
where  6,  is  the  polar  angle  for  the  .yth  mode's  wave  vector 
and  6^  is  a  critical  angle  associated  with  the  cavity  geometry, 
it  is  straightforward  to  show  that 


For  the  numerical  compulations  to  follow,  it  will  be  more 
convenient  to  express  the  length  L  appearing  in  Eqs.  (17)  and 
(20)  in  terms  of  a  zpf-mode  spacing  Ato^.  Thus,  since 
=  277“c*Vw2iAcUj,  these  two  equations  become 


(^?2)+(^n)  =  -^^^^i/^2i|-2  C;(x  i,S-^s 


X  I  (crp(/'))exp[/5,(r-/')]tJ/', 

Jo 

(2  la) 


( d-S-,)  +  ( d-^o'^)  =  -  1 A<-21  I ' 2  C;{x-tsKf<^s 

'  777?- C*  sk 

X  \\an.(t'))cos[S,ir-t^)]dt\  (21b) 

Jo 


B,  Colored  vacuum 

In  order  to  employ  Eqs.  (21),  some  account  must  be  given 
of  the  vacuum  environment  in  which  the  atoms  find  them¬ 
selves.  Specifically,  the  C,  and  (i-£^x)  must  be  evaluated. 
Here,  as  illustrated  in  Fig.  1,  we  restrict  our  attention  to  an 


C 


=  (I-7?)  1  + 


4R 


1/2 


1 


Here,  R  is  the  mirror  reflectivity  and  d  is  the  cavity  length. 
Alternatively,  C,  equals  unity  for  the  modes  with  |cos(^j)| 
<cos(0r).  The  linewidth  of  the  confocal  etalon  (half  width  at 
half  maximum)  for  the  field  evaluated  at  the  center  of  the 
cavity  is 


c 

Aa>c  =  ^sin 


(1-/?) 

l^R 


(23) 


and  its  free  spectral  range  is  just  did. 

Since  the  above  density-matrix  analysis  has  assumed  that 
the  C,  are  unrelated  to  the  density-matrix  elements  (i.e.,  that 
the  atom's  dynamics  do  not  influence  the  zpf  modal  ampli¬ 
tudes  within  the  cavity),  the  present  theory  requires 
Lo)cly^>\.  where  is  the  spontaneous  decay  rate  in  the 
cavity.  In  other  words,  the  present  theory  assumes  that  the 
time  scale  for  the  zpf  to  reach  equilibrium  within  the  cavity 
is  shorter  than  the  time  scale  for  the  atom  to  perturb  the  zpf 
by  its  (real)  radiative  field.  It  is  to  be  noted,  however,  that 
this  is  not  a  fundamental  limitation  to  the  present  semiclas- 
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sical  approach,  as  in  principle  it  should  be  possible  to  relate 
the  to  the  atom's  dynamics. 

To  determine  (x-  we  consider  a  cavity  orientation  as 
illustrated  in  Fig.  1,  with  the  cavity  axis  parallel  to  z  so  that 
the  laser  field  does  not  couple  into  the  etalon.  Then,  defining 

as  the  wave  vector  of  the  5  th  zpf  mode,  we  have 

1^'/  . 

=  C;(  6',)[  1  -  sin-(  e,}cos-(  9  J] , 

(24) 

where  is  the  azimuthal  angle  for  the  5th  mode's  wave 
vector  and  we  have  specifically  indicated  the  dependence 
on  6,. .  Averaging  Eq.  (24)  over  477  sr  [26], 


{x-i,S-=cl{es) 

X 


477  J 

COS-^(^c) 


.  3(k— 1) 

- : - 


COS(  ■ 


(25) 


Thus,  in  the  colored  vacuum  of  a  confocal  etalon  Eqs.  (21) 
become 


+  = - 'irkp~  I  1 1  ’  Jq  ^  ‘  ^ ^  ^ 


dr'  (26a) 


and 


( 4.)  +  (<^^)  =  -  ^22(^' )) 

2  ’7.t‘^iCos[<5^.(r-f')]  dr'.  (26b) 


Averaging  these  over  the  time  step  ST,  we  obtain 


fc(0 


-V 


:%( 0)21+5,)  .  /  5,5r\ 


5,5r 


-sinl 


-  cos(  5,f )  (28a) 


and 


277,(^21+^.)  , 

F,(t)  =  Z  - - sm| 


5,5r 


5,5r\ 


sin(5,?),  (28b) 


which  are  the  field  functions  actually  employed  in  the  nu¬ 
merical  simulations. 

The  effect  of  this  averaging  process  is  to  multiply  each  of 
the  zpf  modes  by  a  sine  function  [i.e.,  sinc(x)  =  sin(x)/x], 
which  is  just  the  Fourier  transform  of  the  discrete  time  step 
employed  in  the  simulation.  We  can  therefore  limit  the  range 
of  modes  in  the  summation  by  choosing  the  two  values  of  S^. 
where  the  sine  function  has  its  first  zero,  since  for  larger 
values  of  \S^\  the  averaging  process  causes  the  modes  to 
have  negligible  amplitude.  Thus,  to  compute  the  zpf  func¬ 
tions  we  need  only  specify  the  modal  spacing  and  the 
time  step  ST,  and  once  computed  the  zpf  functions  are  in¬ 
corporated  into  the  density-matrix  equations  just  like  any 
other  discretized  time  function.  Simulations  are,  of  course, 
repeated  for  various  values  of  A  and  ST  to  ensure  that  the 
results  do  not  depend  on  the  specific  numerical  parameters. 

D.  Full  density-matrix  evolution 

Combining  Eqs.  (2),  (26),  and  (28),  we  obtain  the  semi- 
classical  density-matrix  equations  for  the  evolution  of  a  two- 
level  atom  in  the  presence  of  a  real  (possibly  strong)  field  in 
a  colored  vacuum: 


2a>2iAa>j 

&22=-^1  - |a12iI 


X 


j  (T'y^{t' 

Jo 


)F,(t-r)dr 


(29a) 


and 


C.  Numerical  considerations  for  density-matrix  computations 

To  employ  Eqs.  (26)  in  a  numerical  simulation,  account 
must  be  taken  of  the  discrete  time  step  employed  in  the  simu¬ 
lation  and  the  necessity  for  employing  a  finite  number  of  zpf 
modes.  In  our  paper,  we  first  define  a  time  scale  T  over 
which  we  wish  to  simulate  the  atomic  dynamics,  and  then 
discretize  this  into  elements  of  length  ST.  The  ST  are  small 
enough  so  that  the  density-matrix  elements  do  not  change 
significantly  over  this  interval.  We  then  define  the  two  field 
functions  and  F,(t)  that  are  required  by  Eqs.  (26): 


iFl  (o;,AcUy 

‘^i2=Y^‘M2<r3,-l)--^l)a2.|- 


X 


\‘(TAr')[F,{t-F)^iF,{t-t')]dt'.  (29b) 
Jo 


(Here,  and  in  what  follows,  we  drop  the  angular  brackets 
indicating  an  average  over  the  zpf,  this  now  being  under¬ 
stood.)  To  put  these  into  a  more  recognizable  form,  we  let 
so  that 


Fci‘)  =  'Z  %(w21+^.)cos(5,/)  (27a) 

S 


(722=  — n  Im[cri2]“ 


2W2]A  Wy 

rrfic'^ 


Wzx\- 


and 

=  2  ?7v(ct>'>i+ ^.v)sin(4^)*  (27b) 

.V  and 


(T22{f)F,{t-F)dF 


(30a) 
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(Ti2=  “  /A£0-i2+  ^ 


a>2jAa>,v 

Trhc''* 


terms  in  Eq.  (32).  In  the  limit  of  /— and  Ao)^  >0,  the 
modal  sum  of  the  exponential  terms  in  and  k2  equals  zero, 
so  that  we  have 


Jo 


(30b) 


E.  Lamb  shift 


limALamb=^^ 


a>;a)„ 


(34) 


To  include  mass  renormalization  into  Eq.  (34)  [22],  we 
recognize  that  in  the  limit  of  large-frequency  zpf  modes. 


Though  it  may  not  be  readily  apparent,  a  Lamb-shift  term 
for  the  two-level  atom  is  included  in  Eq.  (30b).  To  demon¬ 
strate  this,  we  need  only  compare  the  first  and  last  terms  on 
the  right-hand  side  of  Eq.  (30b  ),  thereby  recognizing  that  the 
imaginary  part  of  the  last  term  must  act  like  a  zpf  detuning 
term: 


5Lamb=^^p5  ’ 

Subtracting  Eq.  (35)  from  Eq.  (34)  then  yields  an  expression 
for  the  observable  Lamb  shift  of  |m). 


0)2 1 A 

7Thc‘^ 


77,  Wj  sin[<5.f(/-/')] 


di' 


A(Og  ^ 

Ai^ntnh^  V  ,3 Vs\f^mm  ^nn 


Trhc'' 


0)7“  o)„ 


(36) 


^^12^  Limb* 


(31) 


However,  in  order  to  obtain  a  more  realistic  evaluation  of  an 
atom's  Lamb  shift,  specifically  a  cavity-induced  change  in 
the  Lamb  shift  that  may  be  compared  to  experiment,  it  is 
necessary  to  go  beyond  the  two-level  approximation  and 
consider  zpf-induced  virtual  transitions  to  other  atomic 
states. 

Our  starting  point  will  be  Eq.  (32)  from  Ref.  [14],  which 
describes  the  zpf-induced  evolution  of  an  excited-state  am¬ 
plitude,  ,  up  to  second  order  in  perturbation  theor>'.  (The 
b,„  are  expansion  coefficients  of  the  atomic  wave  function  in 
terms  of  its  unperturbed  basis  states,  averaged  over  the  zpf.) 
Specifically, 


+  e 


ih' 

Jo 


(32) 


where  m  and  n  refer  to  unperturbed  energy  eigenstates,  ki 
^  (X),- a)„,„  and  k2^io,- .  In  writing  Eq.  (32),  the  origi¬ 
nal  equation  from  Ref.  [14]  has  been  modified  so  as  to  make 
it  applicable  to  a  colored  vacuum;  it  has  been  averaged  over 
477  sr,  and  rapidly  oscillating  terms  in  the  sum  over  states 
have  been  averaged  to  zero.  Recognizing  that  the  imaginary 
term  of  Eq.  (32)  corresponds  to  the  Lamb  sliift  of  j/tt), 
it  is  Straightforward  to  show  that 


Ai^„b=  Iml  2  ' 


iTrhc 

1 

ik  n 


ikf 


(33) 


Though  Eq.  (36)  is  logarithmically  divergent,  and  would 
require  a  high-frequency  cutoff  for  evaluation,  it  should  be 
recognized  that  our  interest  is  not  in  the  Lamb  shift  itself  but 
in  the  change  in  the  Lamb  shift  caused  by  a  colored  vacuum. 
We,  therefore,  define  the  Lamb  shift  change  as 

c  _  A  cavity  _  A  free  space  (‘yn\ 

^Lamb  ^L^b  “^Lamb  *  ^  ' 

Since  77^  =  273  in  free  space. 


(38) 


Ao)^ 


= 


3  2  Vx  ^  ,  I 


where  the  superscript  on  Si^\y  indicates  that  this  value  is 
specific  to  the  energy  eigenstate  |m).  Since  the  cavity  s  mir¬ 
ror  reflectivity  must  go  to  zero  for  either  veiy  low  or  very 
high  frequency  zpf  modes,  77^  asymptotes  to  2/3  in  the  high- 
frequency  limit  thereby  ensuring  the  convergence  of  ^Lamb* 
In  order  to  evaluate  the  Lamb  shift  numerically,  we  can 
break  the  sum  over  zpf  modes  into  two  contributions,  those 
within  a  cavity-free  spectral  range  of  the  atom  s  resonance 
and  those  outside  this  range.  For  the  zpf  modes  outside  the 
range,  we  convert  the  sum  over  modes  to  an  integral,  thereby 
obtaining 


X 


cu.Ao)^ 


+  I  ^L>( (jj)d (o  I , 


(39) 


where  (Oj  and  cu/,  are  the  low-  and  high-frequency  cutoffs  for 
the  cavity  mirror’s  reflectivity,  respectively;  o)  ,  o)*^,  and 
L(a))  are  given  by 


where  the  prime  on  the  sum  indicates  that  only  nonzero  a)”  =  a>„„, - ,  (40a) 

terms  of  A:  I  and  Jt,  are  included,  as  the  null  values  yield  real  2  ’ 
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CO*  =  CO 


FSR 


nm  9 


(40b) 


Lico)=- 


[77(w)-f]w 


CO--CO,, 


and 


Note  that  is  to  be  chosen  so  that  a;y=r= 


(40c) 


(41) 


m.  VERIFICATION  OF  THE  SEMICLASSICAL 
APPROACH 

As  a  first  check  on  the  semiclassical  description's  validity, 
we  compare  the  predictions  of  Eqs.  (30)  and  (39)  with  the 
experiments  of  Heinzen  and  Feld  [16],  who  passed  a  Ba*^^ 
beam  through  an  etalon  and  examined  weak-field  resonance 
fluorescence  while  the  atoms  were  located  between  the  etal- 
on’s  two  concentric  mirrors.  Heinzen  and  Feld  excited  the 
6.v“  ^Sq-6s6p  transition  at  553.5  nm,  and  recorded  the 
width  and  central  position  of  the  fluorescence.  In  this  way, 
they  were  able  to  obtain  both  the  radiative  decay  rate  of  the 
excited  state,  yrad^  change  in  the  transition's  Lamb 

shift  for  a  quantum  system  perturbed  by  a  colored  vacuum. 
Their  mirror's  separation  was  5  cm:  the  mirror  diameter  was 
1.88  cm,  and  the  mirror  reflectivities  were  low.  Heinzen  and 
Feld's  results  are  reproduced  in  Figs.  2(a)  and  2(b)  as  a  func¬ 
tion  of  cavity  detuning  (triangles),  where  zero  detuning  im¬ 
plies  that  the  cavity  resonance  overlaps  the  atomic  reso¬ 
nance.  The  open  circles  in  Fig.  2(a)  correspond  to  the 
predictions  of  Eqs.  (30)  with  0  =  0,  based  on  the  cavity  ge¬ 
ometry  and  mirror  reflectivity  parameters  provided  by  Hein¬ 
zen  and  Feld  [27].  For  the  theory,  it  was  assumed  that  at  time 
r  =  0  the  atom  was  in  a  superposition  of  the  ground  state  and 
the  excited  state,  and  three  decay  rates  were  computed  cor¬ 
responding  to  the  three  independent  density-matrix  elements. 
In  all  cases,  the  coherence  and  population  decays  were  ex¬ 
ponential,  and  the  computed  decay  rates  of  the  coherence 
were  half  that  of  the  population. 

The  change  in  the  Lamb  shift  was  computed  using  Eq. 
(39),  where  we  set  27rc/a;/  =  900  nm  and  IttcIcji, 
=  350nm.  The  mirror  reflectivity  was  a  constant  value  be¬ 
tween  these  wavelengths  and  zero  outside  this  range.  The 
relevant  transitions  for  the  Lamb  shift  calculation  are  shown 
in  Fig.  3,  and  the  corresponding  dipole  moments  are  col¬ 
lected  in  Table  1;  the  indicated  dipole  moments  come  from 
the  Einstein  A  coefficients  of  Niggli  and  Huber  [28].  Two 
theoretical  curves  are  shown  in  Fig.  2(b).  For  the  solid  curve, 
only  the  553.5  nm  transition  was  included  in  the  Lamb  shift 
calculation,  while  the  dashed  curve  corresponds  to  inclusion 
of  all  transitions  shown  in  Fig.  3.  As  might  have  been  antici¬ 
pated  based  on  an  examination  of  the  dipole  moments,  the 
6,v“  ^Sq-6s6p  transition  has  the  greatest  influence  on 
the  Lamb  shift's  change. 

The  agreement  between  theory  and  experiment  shown  in 
Figs.  2(a)  and  2(b)  provides  strong  evidence  for  the  weak- 
field  validity  of  the  present  semiclassical  description  of  ra¬ 


FIG.  2.  (a)  Spontaneous  decay  rate  in  a  cavity  as  a  function  of 
cavity  tuning.  Diamonds  correspond  to  the  data  of  Heinzen  and 
Feld  [16],  while  circles  correspond  to  the  present  theory,  (b)  Lamb 
shift  change  in  the  cavity  as  a  function  of  cavity  tuning.  Again, 
diamonds  correspo;id  to  the  data  of  Heinzen  and  Feld.  For  the  solid 
line,  only  the  6s'  '5o‘656p  ‘P|  transition  was  included  in  the  Lamb 
shift  computation,  while  in  the  case  of  the  dashed  line  all  transitions 
shown  in  Fig.  3  were  included. 


diative  decay  and  the  Lamb  shift,  especially  when  it  is  rec¬ 
ognized  that  no  free  parameters  were  employed  in  the 
density-matrix  equations.  However,  what  remains  to  be  as¬ 
sessed  is  the  semiclassical  approach's  validity  in  the  case  of 
strong  field/atom  interactions.  In  this  regard,  the  present 
theory  is  most  amenable  to  a  computation  of  transient  nuta¬ 
tion  in  a  colored  vacuum  [29].  Unfortunately,  no  such  experi¬ 
ments  have  been  done,  at  least  in  cavities  conforming  to  the 
present  theory’s  limitations  on  cavity  linewidth.  Conse¬ 
quently,  we  have  used  the  present  theory  to  compute  tran¬ 
sient  nutation  in  free  space,  and  for  added  rigor  have  consid¬ 
ered  nutation  induced  by  a  stochastic  laser  field  (i.e.,  a  phase 
diffusion  field  or  PDF).  We  compare  our  results  to  well- 
established  expectations  [30]. 

For  a  PDF  of  lin width  jp  full  width  at  half  maximum,  it 
is  well  known  that  the  field's  phase  fluctuations  enhance  the 
coherence's  relaxation  rate  while  leaving  the  population  re¬ 
laxation  rate  unaffected.  Specifically,  the  transient-nutation 
signal  for  a  two-level  atom  becomes 


/  e  (f>)\ 

a22(t)  =  1  j , 


(42) 


where  (7o=  corresponds  to  the  steady-state  value  of  the 
excited-state  population.  Also, 


013815-7 


JAMES  CAMPARO 


PHYSICAL  REVIEW  A  65  013815 


A-^i-l-  y^ 

'2  ’ 

(43a) 

(43b) 

and 

</>  =  tan“*  ^  j .  ^43c) 

Here,  corresponds  to  the  coherence  decay  rate,  which 
according  to  well-established  theory  and  experiment  should 
equal  (A  21+  rr)/2,  and  A  21  is  the  spontaneous  decay  rate  in 
free  space. 

In  the  test,  we  again  considered  laser  excitation  of  the 
6,v"  ^So-6s6p  ^Pi  barium  transition.  The  PDF  was  numeri¬ 
cally  simulated  for  some  given  value  of  and  the  density- 
matrix  equations  were  then  solved  numerically  using  a  vari¬ 
able  step-size  Runge-Kutta-Fehiberg  method  [31].  The 
computed  transient- nutation  signal  was  averaged,  and  then 
using  a  nonlinear  least  squares  procedure  [32]  was  fit  to  Eq. 

TABLE  I.  Selected  transition  dipole  moments  for  barium.  Val¬ 
ues  are  inferred  from  the  lifetime  measurements  of  Niggh  and  Hu¬ 
ber  [28]. 

Transition  Dipole  moment,  \p\  (esu  cm) 


6.v6/j  '/’|-6.v"  'Sq 

8.0X  10' 

6.v8.v  '5(i-6.v6/>  'Pi 

2.1X  10' 

(>sld  'D:-6.v6/7 

3.6  X  10' 

6p- ^PQ-6.t6p  'P, 

2.8  X  10 

5d6p  '5o 

5.8  X  10 

FIG.  4.  Atomic  dephasing  rate  in  free  space  as  a  function  of 
phase-diffusion-field  (PDF)  linewidth.  The  solid  line  corresponds  to 
expectations  based  on  experiment  and  theory.  The  open  circles  cor¬ 
respond  to  the  results  of  the  present  theory. 

(42)  with  yo  as  the  one  free  parameter.  The  results  of  the  test 
are  shown  in  Fig,  4  as  open  circles,  where  the  value  of  y2 
from  the  numerical  simulation  is  plotted  as  a  function  of  the 
input  value  of  yp .  The  solid  curve  in  the  figure  corresponds 
to  the  expected  value  of  y2.  (A21+ yr)/--  The  excellent 
agreement  between  the  numerically  simulated  value  of  y2 
and  (A2i+y/r)/2  not  only  lends  confidence  to  the  present 
theory's  ability  to  simulate  strong  field/atom  interactions,  it 
also  suggests  that  the  semiclassical  description  of  radiative 
decay  may  find  application  in  simulating  the  effects  of  col¬ 
ored  vacua  in  studies  of  the  stochastic-field/atom  interaction. 

IV.  TRANSIENT  NUTATION  IN  A  COLORED  VACUUM 

Given  the  previous  section's  evidence  for  the  validity  of 
the  present  semiclassical  description  of  radiative  decay,  we 
now  consider  the  decay  rate  of  transient  nutation  in  a  colored 
vacuum.  For  the  present  study,  we  confine  our  attention  to  a 
monochromatic  laser  tuned  to  barium  s  6^“  ^Sq-6s6p  P[ 
transition,  and  the  cavity  geometry  of  Heinzen  and  Feld  [16]. 
As  discussed  above,  when  the  cavity  detuning  from  the 
atomic  resonance  A^av^  equals  zero,  the  population  and  co¬ 
herence  decay  rates  in  the  weak-field  limit  are  at  a  maxi¬ 
mum,  with  a  value  determined  by  the  cavity  s  mirror  reflec¬ 
tivity.  Consequently,  based  on  Eq.  (43a)  it  was  originally 
thought  that  the  transient-nutation  decay  rate  would  be  maxi¬ 
mized  for  Acav=0  with  a  value  of  3  y^J^.  Using  the  present 
work's  weak-field  results  for  y^ad .  this  prediction  is  shown  as 
the  solid  line  in  Fig.  5  as  a  function  of  mirror  reflectivity. 

Transient-nutation  decay  in  a  cavity  was  investigated  by 
numerically  solving  Eqs.  (30a)  and  (30b)  for  Aj.av=0  using 
the  parameter  values  listed  in  Table  11.  In  order  to  determine 
the  decay  rate  F,  we  Fourier  transformed  the  transient  nuta¬ 
tion  signal  and  determined  the  width  of  the  Fourier  spectral 
component  at  the  nutational  frequency.  In  this  way  the  deter¬ 
mination  of  r  was  independent  of  any  assumption  regarding 
the  functional  form  of  the  nutation  signal  in  the  cavity,  as 
would  have  been  the  case  had  we  fit  the  nutation  signal  to 
Eq.  (42).  The  transient-nutation  decay  rates  obtained  in  this 
way  are  shown  as  circles  in  Fig.  5.  What  is  strikingly  appar¬ 
ent  is  that  the  numerical  results  indicate  a  decay  rate  signifi¬ 
cantly  different  from  the  expected  3  values,  and  that 
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FIG.  5.  Tran.sienHiutation  decay  rate  F  in  a  cavity  as  a  function 
of  the  cavity  mirrors'  reflectivity.  The  solid  line  corresponds  to  an 
anticipated  value  based  on  the  extrapolation  of  weak- field  cavity 
decay  rates  into  the  strong-field  regime.  Circles  correspond  to  the 
results  of  the  pre.seni  theory  for  Acav“0'  while  diamonds  corre¬ 
spond  to  the  present  theory  for 


atom  appears  to  “filter''  the  vacuum  field  fluctuations,  pro¬ 
ducing  an  output  that  forces  the  excited  state  to  decay. 

Taking  advantage  of  the  convolution  theorem  [33],  Eq. 
(44)  can  be  reexpressed  as 

tt>2iA(u_y  ^ 

cr 22^  ^  Im[  0^12]  ^ 

f  +» 

X  Z22{<o)^cMe‘'"'dw.  (45) 

J  -CO 

where  and  4E>c  are  the  Fourier  transforms  of  0-22  and  F^. , 
respectively.  Given  the  oscillatory  form  of  it  is  clear 
that  1 22  is  a  sum  of  S  functions: 

2^22(^)  “ (46) 

which  on  substitution  into  Eq.  (45)  and  recognizing  that 
tl>,.(cu)  =  <l>^(- w)  yields 


for  mirror  reflectivities  less  than  0.5  the  nutational  decay  rate 
is  actually  less  than  its  free-space  value.  If  the  cavity  is  de¬ 
tuned  from  the  atomic  resonance  by  the  Rabi  frequency,  then 
the  decay  rates  shown  as  diamonds  in  Fig.  5  are  obtained.  In 
this  case,  transient-nutation  decay  in  the  cavity  is  faster  than 

To  gain  a  semiquantitative  understanding  of  the  nutational 
decay  rate's  sensitivity  to  cavity  detuning  (and  also  to  high¬ 
light  the  semiclassical  approach's  intuitive  value),  we  rewrite 
Eq.  (30a)  for  t>  1/A<u^. ,  so  that  the  upper  limit  of  integration 
may  be  replaced  by  4^.  Moreover,  since  we  are  considering 
nutation,  to  first  approximation  we  can  assume  that  the 
strong  field  causes  0*22  to  oscillate  at  the  Rabi  frequency  [i.e., 
cr22=cos(no].  This  then  yields 


<T22—  Im[(7i2]“" 


0'22i^'~  T)Fc(T)dT^ 
(44) 


where  the  integration  variable  has  been  redefined.  In  this 
form,  the  second  term  on  the  right-hand  side  of  Eq.  (44)  has 
the  appearance  of  a  linear  system's  output  [33];  that  is,  the 


TABLE  II.  Parameters  used  to  compute  transient  nutation  in  a  . 
colored  vacuum. 


Parameter 

Value 

u’:2(^  =  0),  cr,2(r=0) 

0 

Resonant  wavelength 

553.5  nm 

Excited-state  decay  rate  in  free  space,  A  21 

1.19X  la'^sec-' 

Transition  dipole  moment,  I/X21I 

8.02X  10‘'*esu  cm 

Cavity  length 

5  cm 

Mirror  diameter 

1.88  cm 

Laser  intensity 

100  W/cm- 

Rabi  frequency,  fl 

6.97X10^sec“' 

Laser  detuning  from  atomic  resonance.  A^ 

0 

Simulation  time  scale,  T 

100/A  21 

zpf  time-scale  discretization,  ST 

r/8192 

Im[cri2]“rcr22(0,  (47) 


where 


(x)l,Ao)s , 

(48) 

(Note  that  in  the  limit  of  zero  cavity-mirror  reflectivity,  F 
^4w'^^\jLi2[\‘/3^ic^^  the  Einstein  A  coefficient  for  |2)  to  |l) 
spontaneous  decay  in  free  space.) 

Thus  in  the  semiclassical  approach,  the  nutational  decay 
rate  in  a  colored  vacuum  is  seen  to  arise  from  an  “atomic 
filtering"  of  the  vacuum  field  fluctuations,  with  the  atom  re¬ 
sponding  principally  to  that  Fourier  frequency  component  of 
the  vacuum's  fluctuations  corresponding  to  the  Rabi  fre¬ 
quency.  When  the  cavity  is  tuned  to  the  atomic  resonance, 

0^(0)  and  hence  f  can  be  relatively  small.  However,  when 
the  cavity  is  detuned  from  resonance  by  the  Rabi  frequency, 
th,.(n)  is  at  its  maximum  value  and  nutational  decay  is 
rapid.  It  is  to  be  noted  that  this  conclusion  is  consistent  with 
the  works  of  Keitel  et  al.  [17]  and  Lewenstein,  Mossberg, 
and  Glauber  [34],  who  employed  QED  to  investigate  the 
fluorescence  triplet  in  a  colored  vacuum.  Specifically,  Keitel 
et  al,  found  that  the  linewidth  of  the  triplet's  central  compo¬ 
nent  only  depended  on  the  vacuum  spectral  density  at  the 
Rabi  sidebands. 


V.  SUMMARY 

Given  the  history  of  past  attempts  to  describe  spontaneous 
emission  semiclassically,  I  wish  to  emphasize  that  the  pur¬ 
pose  of  the  present  work  is  not  to  propound  some  rival  to 
QED.  It  has  been  known  for  decades  that  if  matter  is  quan¬ 
tized,  then  consistency  dictates  quantization  of  the  electro¬ 
magnetic  field  [35].  Rather,  the  puipose  of  the  present  work 
is  simply  to  extend  the  range  of  problems  that  can  be  ad¬ 
dressed  tlirough  semiclassical  methods.  Here,  a  semiclassical 
description  of  radiative  decay  has  been  incorporated  into  the 
density-matrix  equations  in  order  to  develop  a  completely 
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semiclassical  methodology  for  dealing  with  strong- field/atom 
interactions  in  a  colored  vacuum. 

One  application  of  this  methodology,  which  could  have 
computational  advantages  over  QED,  would  be  investiga¬ 
tions  into  multiphoton  processes  taking  place  in  a  low-2 
cavity,  possibly  driven  by  a  stochastic  field  [36].  In  this  case, 
the  advantage  pertains  to  the  complete  problem,  and  not  sim¬ 
ply  to  a  calculation  of  vacuum  effects.  For  example,  since  the 
stochastic  character  of  a  laser  is  (arguably)  most  easily  mod¬ 
eled  by  treating  the  laser  as  a  classical  field,  and  thereby 
most  easily  connected  to  laboratory  measurements,  consis¬ 
tency  would  dictate  that  the  vacuum  field  should  also  be 
treated  classically  in  a  stochastic-field/atom  interaction  prob¬ 
lem.  Treating  the  vacuum  field  by  using  QED,  while  accu¬ 
rate,  necessitates  treatment  of  the  stochastic  laser  as  a  quan¬ 
tized  field,  which  is  not  always  the  simplest  approach  for  this 
type  of  problem.  In  fact,  it  was  the  consideration  of  such 
problems  that  prompted  our  interest  in  the  present  work. 

Of  course,  it  must  be  remembered  that  the  present  theory 
is  only  an  approximation  to  QED,  and  therefore  there  are 
bounds  on  its  range  of  validity.  For  example,  by  limiting  our 
discussion  to  the  two-level  atom  we  have  ignored  much  of 
the  vector  nature  of  spontaneous  emission.  While  this  could 
probably  be  accommodated  within  the  semiclassical  theory, 
there  seems  little  to  be  gained  at  the  present  time  by  the 
added  effort.  Of  greater  significance  is  the  present  theory's 
focus  on  spontaneous  decay  and  the  assumption  that  the 
atom  does  not  alter  the  zero-point  field's  spectral  density. 
Specifically,  taking  a  conservative  point  of  view,  we  can  only 
expect  the  semiclassical  theory  to  agree  with  QED  up  to 
second  order  in  the  vacuum's  perturbation.  This  in  turn  im¬ 
plies  that  the  present  theory  should  only  be  trusted  to  de¬ 
scribe  effects  associated  with  second-order  vacuum  field  sta¬ 
tistics.  Finally,  it  must  be  noted  that  in  order  to  treat 
radiation-reaction  in  a  simple  way,  the  density-matrix  equa¬ 
tions  were  averaged  over  the  vacuum  field.  As  a  conse¬ 
quence,  the  density-matrix  treatment  presented  here  is  lim¬ 
ited  to  problems  that  deal  with  one-time  averages  of  the  zpf, 
and  therefore  cannot  account  for  processes  that  are  associ¬ 
ated  with  two-time  averages,  for  example,  the  resonance 
fluorescence  triplet  [37]. 

Before  concluding,  it  seems  worthwhile  to  suggest  that 
the  semiclassical  methodology  discussed  here  might  find  ap¬ 
plication  beyond  atomic  and  optical  physics.  For  example, 
the  present  approach  may  prove  useful  in  studying  open 
quantum  systems  and  the  issues  surrounding  quantum  deco¬ 
herence  [38].  Complementing  quantum-field  treatments  of 
this  problem,  the  present  semiclassical  approach  might  pro¬ 
vide  additional  insights  given  its  different  intuitive  appeal. 
Additionally,  one  could  imagine  using  the  present  approach 
as  a  template  for  developing  a  semiclassical  approximation 
to  quantum  gravity  in  order  to  study  processes  associated 
with  the  gravitational  vacuum.  In  this  vein,  it  is  worth  citing 
the  work  of  Ross  and  Moreau  [39],  who  have  discussed  the 
characteristics  of  a  classical,  gravitational  zero-point  field  in 
the  context  of  an  approximation  to  the  real  gravitational 
vacuum.  Though  both  these  possible  applications  of  the 
present  work  are  highly  speculative,  they  do  highlight  the 


benefits  to  be  gained  by  extending  semiclassical  methodolo¬ 
gies  as  far  as  possible. 

acknowledgment 

This  work  was  supported  under  U.S.  Air  Force  Contract 
No.  F04701-93-C-0094. 


APPENDIX 

In  Ref.  [14],  the  Schrodinger  equation  for  an  atom  inter¬ 
acting  with  the  zero-point  field  (zpf)  of  Eq.  (3)  and  its  own 
radiation-reaction  field  [i.e.,  Eq.  (7)]  was  solved  up  to  second 
order  in  perturbation  theory.  In  typical  fashion,  the  atomic 
wave  function  \cp)  was  expanded  in  terms  of  eigenfunctions 
of  the  unperturbed  Hamiltonian  |/?)  and  the  evolution  of  the 
expansion  coefficients  (or  wave-function  amplitudes)  a,j{t)^ 
under  the  influence  of  the  two  field-atom  interactions  was 
derived.  Specifically,  using  overbars  in  tliis  appendix  to  in¬ 
dicate  an  average  over  the  zpf,  it  was  found  that 


|9)  =  X  (Al) 

n 

with  dfj  the  sum  of  a  zpf  term,  and  a  radiation-reaction 
term  a^^  such  that 


X?  1  12  3 


(A2a) 


IA„/X 


a>; 


(A2b) 


Of  significance  is  the  fact  that  the  two  rate  terms  interfere 
constructively  to  produce  spontaneous  decay,  and  interfere 
destructively  to  inliibit  spontaneous  excitation. 

In  order  to  create  a  density-matrix  that  includes  this 
constructive/destructive  interference  aspect  of  the  atom  s  in¬ 
teraction  with  the  zpf  and  its  own  radiation-reaction  field,  we 
would  need  to  employ  the  zpf-averaged  wave  function  of  Eq. 
(Al): 


=  X  ''“"’"IwXnl.  (A3) 

nw 

Then,  for  a  two-level  atom  it  is  straightforward  to  show  that 
in  a  rotating  frame,  where  pij—aijc 

r 

r<722  — 

^  r 

-  7  ^^21  ^  ^22 

Here,  F  is  the  spontaneous  decay  rate  (i.e., 
4l^^l|'^02,/3/ic^),  and  we  have  normalized  the  density  ma¬ 
trix.  (Since  the  radiation -reaction  perturbation  of  Ref.  [14]  is 
non-Hermitian,  population  from  the  excited  state  is  lost  un¬ 
less  the  requirement  of  normalization  is  imposed  on  the  den¬ 
sity  matrix.)  As  is  obvious  from  Eq.  (A4),  construction  of  the 
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density  matrix  from  the  zpf-averaged  wave  function  yields 
the  proper  spontaneous  decay  behavior. 

In  a  density-matrix  treatment  of  stochastic  processes, 
however,  it  is  customary  to  employ  the  atomic  wave  func¬ 
tions  prior  to  stochastic  averaging.  The  ensemble  average 
behavior  of  the  atomic  system  is  then  obtained  by  averaging 
the  density-matrix  elements  over  the  stochastic  process  [40]. 
In  this  fashion,  which  is  essentially  that  of  the  present  work, 
we  obtain 


=  2  a  mate  (A5) 

mn 

or  in  the  rotating  frame 


cr<->s|<p)(^l  =  2  a„at\m){n\.  (A6) 

mn 

The  problem  with  for  describing  spontaneous  emission 
semiciassically  is  that  the  constructive/destructive  interfer¬ 
ence  aspects  of  the  zpf  and  radiation-reaction  perturbations 
will  not  be  manifested  in  the  terms  if  the  perturbations 
of  Ref.  [14]  are  employed.  Basically,  if  we  are  to  inhibit 
spontaneous  excitation  through  the  interference  mechanism, 
then  we  need  to  average  over  the  wave-function  amplitudes 
and  not  the  probabilities  for  an  atom  being  in  a  specific  state. 

An  expedient  solution  to  this  density-matrix  averaging 
ambiguity  lies  in  the  fact  that  the  radiation-reaction  pertur¬ 
bation  is  not  gauge  invariant,  which  implies  that  it  is  only  a 
computational  tool.  Therefore,  we  have  some  liberty  to  de¬ 
fine  this  perturbation  so  as  to  make  the  evolution  of 
consistent  with  the  evolution  of  Considering  for  the 


moment  just  the  diagonal  components  of  the  density  matrix 
we  can  write  their  evolution  as  the  sum  of  a  zpf- 
induced  evolution  and  a  radiation-reaction  evolution: 


(A7) 


We  then  employ  our  freedom  on  the  form  of  the  radiation- 
reaction  perturbation  to  write  and  d^^a'o 

=  which  is  just  the  ansatz  of  Eq.  (9).  In  this  way, 

spontaneous  excitation  of  the  ground  state  is  inhibited  while 
spontaneous  decay  of  the  excited  state  proceeds  with  the 
proper  rate. 

In  order  to  show  that  the  above  choice  of  radiation- 
reaction  evolution  leads  back  to  Eq.  (A4),  we  can  consider 
Eqs.  (26)  in  the  case  of  free  space,  where  77=2/3.  In  that 
situation,  since  the  sum  over  zpf  modes  will  be  dominated  by 
terms  with  =  we  can  let  07^^0721  so  that 


277  Jq  j  =  0 

r  °° 

- f  cos[<5^(r-;')]Aw^.dr'. 


(A8b) 


Changing  the  sum  over  modes  to  an  integral,  and  ignoring 
terms  associated  with  the  Lamb  shift,  Eqs.  (A8)  yield  &\2 
=  and  0*22^=  “rcr22\  which  with  normalization 

yields  a  density-matrix  evolution  that  is  equivalent  to  the 
^  of  Eq.  (A4). 
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LABORATORY  OPERATIONS 


The  Aerospace  Corporation  functions  as  an  “architect-engineer”  for  national  security  programs,  spe¬ 
cializing  in  advanced  military  space  systems.  The  Corporation’s  Laboratory  Operations  supports  the 
effective  and  timely  development  and  operation  of  national  security  systems  through  scientific  research 
and  the  application  of  advanced  technology.  Vital  to  the  success  of  the  Corporation  is  the  technical 
staffs  wide-ranging  expertise  and  its  ability  to  stay  abreast  of  new  technological  developments  and 
program  support  issues  associated  with  rapidly  evolving  space  systems.  Contributing  capabilities  are 
provided  by  these  individual  organizations: 

Electronics  and  Photonics  Laboratory:  Microelectronics,  VLSI  reliability,  failure  analy¬ 
sis,  solid-state  device  physics,  compound  semiconductors,  radiation  effects,  infrared  and 
CCD  detector  devices,  data  storage  and  display  technologies;  lasers  and  electro-optics,  solid 
state  laser  design,  micro-optics,  optical  communications,  and  fiber  optic  sensors;  atomic 
frequency  standards,  applied  laser  spectroscopy,  laser  chemistry,  atmospheric  propagation 
and  beam  control,  LIDAR/LADAR  remote  sensing;  solar  cell  and  array  testing  and  evalua¬ 
tion,  battery  electrochemistry,  battery  testing  and  evaluation. 

Space  Materials  Laboratory:  Evaluation  and  characterizations  of  new  materials  and 
processing  techniques:  metals,  alloys,  ceramics,  polymers,  thin  films,  and  composites; 
development  of  advanced  deposition  processes;  nondestructive  evaluation,  component  fail¬ 
ure  analysis  and  reliability;  structural  mechanics,  fracture  mechanics,  and  stress  corrosion; 
analysis  and  evaluation  of  materials  at  cryogenic  and  elevated  temperatures;  launch  vehicle 
fluid  mechanics,  heat  transfer  and  flight  dynamics;  aerothermodynamics;  chemical  and 
electric  propulsion;  environmental  chemistry;  combustion  processes;  space  environment 
effects  on  materials,  hardening  and  vulnerability  assessment;  contamination,  thermal  and 
structural  control;  lubrication  and  surface  phenomena. 

Space  Science  Applications  Laboratory:  Magnetospheric,  auroral  and  cosmic  ray  phys¬ 
ics,  wave-particle  interactions,  magnetospheric  plasma  waves;  atmospheric  and  ionospheric 
physics,  density  and  composition  of  the  upper  atmosphere,  remote  sensing  using  atmos¬ 
pheric  radiation;  solar  physics,  infrared  astronomy,  infrared  signature  analysis;  infrared  sur¬ 
veillance,  imaging,  remote  sensing,  and  hyperspectral  imaging;  effects  of  solar  activity, 
magnetic  storms  and  nuclear  explosions  on  the  Earth's  atmosphere,  ionosphere  and  magne¬ 
tosphere;  effects  of  electromagnetic  and  particulate  radiations  on  space  systems;  space 
instrumentation,  design  fabrication  and  test;  environmental  chemistry,  trace  detection; 
atmospheric  chemical  reactions,  atmospheric  optics,  light  scattering,  state-specific  chemical 
reactions  and  radiative  signatures  of  missile  plumes. 

Center  for  Micro  technology:  Microelectromechanical  systems  (MEMS)  for  space 
applications;  assessment  of  microtechnology  space  applications;  laser  micromachining; 
laser-surface  physical  and  chemical  interactions;  micropropulsion;  micro-  and  nanosatel¬ 
lite  mission  analysis;  intelligent  microinstruments  for  monitoring  space  and  launch  sys¬ 
tem  environments. 

Office  of  Spectral  Applications:  Multispectral  and  hyperspectral  sensor  development; 
data  analysis  and  algorithm  development;  applications  of  multispectral  and  hyperspectral 
imagery  to  defense,  civil  space,  commercial,  and  environmental  missions. 


